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How does a rapidly rotating condensed Bose gas carry extreme amounts of 
angular momentum? The energetically favored state of a not-too-rapidly 
rotating Bose condensed gas is, as observed, a triangular lattice of singly 
quantized vortices. This paper describes the fates of the vortex lattice in 
both harmonic and anharmonic traps when condensates are rotated extremely 
rapidly. 
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1. INTRODUCTION 

Superfluids respond to rapid rotation by forming triangular arrays of 
singly quantized vortex lines. Compared with superfluid '^He,- the proper- 
ties of vortices in Bose-Einstein condensed atomic gases have proven to be 
considerably more accessible. Following the creation of single vortex lines 
in atomic condensates by phase imprinting^ large arrays have been created 
via mechanical rotation, by stirring the condensate with a laser beamj^ 
and by particle evaporation.!^ Upwards of 300 vortex lines have been created 

in rapidly rotating condensates, and properties such as the modes of the 
lattic(P3 

and vortex core structur^^ have been studied in detail. 
The superfluid velocity, v , obeys the quantization condition, § v ■ di = 
{h/m)N^{C), where the line integral is along a closed contour surrounding 
N^{C) singly quantized vortices, h is Planck's constant, and m the particle 
mass. In the neighborhood of a single line, the velocity is in the azimuthal 
direction and has magnitude h/mp, where p is the distance from the line. 
A system containing many vortex lines appears to rotate uniformly with an 
average angular velocity 0,, which is simply related to the (two-dimensional) 
density of vortex lines, n„, by the quantization condition: ^1 = nhn^/m. As 
the system rotates more and more rapidly, the lines become closer and closer. 
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The question I would like to discuss here is what eventually happens to the 
vortex array. In other words, how does a rapidly rotating Bose condensate 
carry large quantities of angular momentum per particle? 

One may be tempted to argue by analogy with a Type II superconductor 
which in the presence of a magnetic field above a critical value, Hd, contains 
an array of vortex lines with quantized flux. As the field is increased the 
line density grows until the cores begin to overlap, at a critical field, Hc2, at 
which point the system turns normal. However, a low temperature rotating 
bosonic system does not have a normal phase to which it can return.^ Un- 
like in a weakly interacting superconductor, where condensation occurs as 
a small dynamical decoration on top of a normal state, Bose condensation 
occurs kinematically. A Bose system must respond differently than a Type 
II superconductor. In He II, the question has never been of experimental 
interest, since the intervortex spacings are typically macroscopic, while the 
core sizes are of order Angstroms, and thus the critical rotation rate, flc2, 
at which the vortex cores would approach each other is unobservably large, 
~ 10^^ rad/sec.^In a weakly interacting atomic condensate, the radius of 
a single vortex core is of order = l/V^irna, where a is the s-wave scat- 
tering length, and n is the particle density. Typically, ~ 0.2 fim, so that 
cores would begin to touch at Qci ~ Sn'^nah/m ~ 10'^ — 10^ rad/sec, an 
experimentally accessible rate. 

The ground state of the system rotating at angular frequency Q, about 
the z axis is determined by minimizing the energy in the rotating frame, 
E' = E — rJLx, where is the component of the angular momentum of the 
system along the rotation axis: 



E' = / d^r 



2m 



1 n2„2 \ L;,|2 1 



+ ( V{r) - -mn^ij 1^1^ + -g[ 



(1) 



where tp is the condensate wave function, ri = {x,y), V{f) is the trapping 
potential, and g = Airh'^a/m determines the strength of the interparticle 
interaction. 

The fate of a rapidly rotating Bose system depends on how the system 
is confined. In typical condensate experiments the system, rotating about 
the z axis, is confined in a harmonic trap of the form, 

V{r^, z) = Kn (uolrl + uj,z^) . (2) 

In this case the centrifugal potential, — ^mO^r^, tends to cancel the trans- 
verse trapping potential, and the system cannot rotate faster than lv± with- 
out becoming untrapped. As — > uj±, the system flattens out, becomes 
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almost two dimensional, and eventually enters quantum Hall-like states. On 
the other hand, in the presence of an anharmonic transverse trapping po- 
tential that grows faster than quadratic, e.g., 

^±(r±) = ^m^iri(l + Ari), (3) 

the system can rotate faster than since it is contained by the anharmonic 
part of the potential. 



2. THE VORTEX CORES SHRINK WITH INCREASING ft 

In fact there is never a phase transition associated with the vortex cores 
overlapping in a rotating Bose condensate. Rather, as the intervortex spac- 
ing becomes comparable to the mean field coherence length, the vortex 
cores begin to shrink, and eventually the core radius scales down with the 
intervortex spacing. 
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Fig. 1. (a) Mean vortex core area as a fraction of the area per vortex vs. the 
rotational rapidity (see text), (b) Measured vortex core area as a fraction of 
the area per vortex vs. 20. /gn. Adapted from Ref. 



In Ref. .11 , and more fully in Ref. 1121 we approached the problem by 
treating the core radius, ^, as a variational parameter. By integrating out 
the short range structure on scales of the vortex separation, one derives an 
effective Hamiltonian, 



E' 



J d^rn{f) 



^2 



! I 

-mr^ + V{r) + afiQ. + -bgn{f 



(4) 
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where n{r) is the smoothed density; the parameter a(^) takes into account 
the local vortex kinetic energy, including the curvature of the order pa- 
rameter, and describes the renormalization the coupling constant, g, 
due to density fluctuations about the smoothed density within each vortex 
cell. The core size is determined by minimizing Eq. Q with respect to ^. 
Figure la shows the resultant mean square area of the vortex core, in a 
harmonic trap, measured in units of the area per vortex line, for the quali- 
tatively accurate model in which the condensate wave function has the form 
ip{p) ~ p, for p < and constant for ^ < p < where I is the radius of 
the (cylindrical) Wigner-Seitz cell around a given vortex: = 1/mQ.. In 
this model, the mean square core area divided by the area per vortex, A, is 
^■^/3£^. The horizontal axis in Fig. 1 is the rotational rapidity,^?/, defined 
by ^/lo±_ = tanhy, a convenient variable for spreading out the region where 
Q.^u}±. The linear rise at small occurs because the core size remains 
constant, while l"^ decreases linearly with The flattening of A with 

increasing is a consequence of the vortex radius scaling with the inter- 
vortex spacing. The upper line shows the exact limit as O ^ lo^. Recent 
JILA measurement^^ of A (as a function of 2h^}/gn), Fig. lb, nicely show 
the expected initial linear rise, followed by the predicted scaling of the core 
radius with intervortex spacing. 



3. LOWEST LANDAU LEVEL REGIME 

As the rotation rate in a harmonic trap approaches uj±, the centrifugal 
potential basically cancels the transverse trapping potential; the cloud flat- 
tens out, and becomes a effectively two dimensional system. Because the 
density n of the system drops, the interaction terms ~ gn, become small 
compared with HQ, Then, as one sees from Eq. with the dynamics is 
that of a particle feeling the Coriolis force alone, a system formally analogous 
to a particle in a magnetic field. HcP predicted that in this limit particles 
should condense into the lowest Landau level (LLL) in the Coriolis force, 
similar to charged particles in the quantum Hall regime. When 2gn <^ Q, 
the states in the next higher Landau level are separated by a gap ~ 2uj±. 
This insight has led to extensive experimental studies in which rotation rates 
in excess of 0.99 ui± have been achieved.^ 

The single particle wave functions in the lowest Landau level have the 
form, (t>fj.{f±) ~ (^/^e~''j-/^'^J-, where ( = x + iy, p = 0,1,2,..., and the 
transverse oscillator length, is given by y^h/mujj_. The LLL condensate 
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wave function is a linear superposition of such states: 

N 

<^LLL(r-l) ~ S^c^C'^e-'^'/24 ~ - Ci) e-^'/2'^i, (5) 

1=1 

where the polynomial Ti^c^C,^ is written as a product over its zeroes, Cij 
which are simply the positions of the vortices in the condensate. The state 
in the regime gn <^ HQ is a direct continuation of the state in the slowly 
rotating regime, HQ <^ gn. 



4. DENSITY PROFILE AND LATTICE DISTORTION 

As long as the total number of vortices is much larger than unity, the 
energy of the cloud in the LLL limit is given 

E' = QN + J d^r{{uj^ - fi)^n(r ) + yn(f (6) 

plus terms involving the trapping potential in the z direction. Here n(r ) 
is the smoothed density profile, = (|(/>lllP); the brackets denote the long 
wavelength smoothing. The energy Q is minimized when the cloud assumes 
a density profile of the Thomas-Fermi form, - an 

inverted parabola - where R is the transverse radius of the cloud. For 
Na/dz S> 1, where is the axial oscillator length, the structure in the 
radial direction will be Thomas-Fermi at large fi, even if it is Gaussian at 
small ri.^ In experiment the density profile indeed remains an inverted 
parabola as Q ^ u;±. 

Since the energy © depends only on the smoothed density, the vortices 
must adjust their locations in order that the smoothed density be an inverted 
parabola. From the arguments in Ref. 1131 we find the relation between the 
smoothed density and the mean vortex density, n^{r±), 

1 1 

-V^lnn(r_L) = + 7rn(r_L). (7) 

For a Gaussian density profile, the vortex density is constant. However, for 
a Thomas-Fermi profile. 

Since the second term is of order compared with the first, the density 

of the vortex lattice is basically uniform (and the vortex array forms an 
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almost perfect triangular lattice). Turning the argument around, very small 
distortions of the vortex lattice from triangular can result in large changes 
in the density distribution Recent measurements of the (percent scale) 
distortions of the vortex lattice at relatively low rotation ratesPare in good 
agreement with theory. 



5. LATTICE MODES 

The collective modes of the vortex lattice in a rotating superfluid have 
been of considerable interest since the 1960's, when Tkachenko^^ showed 
that the lattice supports an elliptically polarized mode, a mode observed in 
helium in 1982.^ To describe the Tkachenko modes, Ref. 1201 reformulated 
the hydrodynamics of rotating superfluids to take into account the elasticity 
of the vortex lattice (including normal fluid, dissipation, and line bending - 
Kelvin - oscillations of the vortex lines in three dimensions); effects of the 
oscillations of the vortex lines at finite temperature on the long range phase 
correlations of the superfluid were discussed in Ref. 1211 In superfluid helium, 
where is always much smaller than characteristic phonon frequencies, sk ~ 
s/R, where s is the sound velocity and R the transverse size, the modes 
frequencies are linear in the wave vector, ujt = {2C2/mn)^^'^k, where C2, the 
shear modulus, equals hn^l/8. However, in Bose condensates it is possible 
to go to the rapidly rotating regime^ ^ ^ where Q > sk. For general Q, the 
in-plane Tkachenko mode frequencies are given by,'^^'^ 

2 ^ 2^ 

'^^ nm {402 + (52 + + C2)/nm)k-^} ' 

where Ci is the compressional modulus of the vortex lattice. Note that at 
low k, the mode frequency becomes ~ fc^, addition the system supports 
a gapped sound mode of frequency, 

a;f = in^+(s^ + —{Ci + C2)) k\ (10) 
V nm / 

Figure 2a shows the frequencies as a function of k, while Fig. 2b shows the 
lowest mode frequency (of wave vectoi^ 5.45/i2) evaluated for a trapped 
rotating ^^Rb gas of = 2.5 x 10^ particles.!^ In experiment,'^ the particle 
number varies from run to run; here we have multiplied the experimental 
frequencies by a factor {N/2.5 x 10^)"^^^ to compare with theory, calculated 
for N = 2.5 X 10^. As one approaches the LLL regime, the shear modulus 
decreases t(P^C2 ~ (81/807r4)ms2 n; the Tkachenko modes are then softer 
by a factor (9/7T^){ms^ /Wn)^/^; see solid curve in Fig. 2a. This softening 
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Fig. 2. a) Gapped sound and Tkachenko mode frequencies vs. wavevector, 
for <C ms^. (b) Frequency of the lowest Tkachenko mode, with the mean 
sound velocity, a decreasing shear modulus, C2 (solid curve) and a constant 
C2 (dashed curve). The data (triangles) are from Ref. El The inset shows 
u)T over the entire range of fi, at constant C2. 



of the modes has recently been measured,'^^ and is in good agreement with 
theory. 

The softness of the Tkachenko modes in the rapidly rotating regime 
leads to infrared singular behavior in the vortex transverse displacement- 
displacement correlations at finite temperature, and in the order parameter 
phase correlations even at zero temperature.l^^SlESl xhe zero point oscil- 
lations of the vortices cause quantum melting of the lattice when Ny ~ 
10 — 20A^, where A'^ is the total number of vortices present.'^^'^ In a finite 
system the single particle density matrix, {ip{r)ip^ {r')), falls algebraically as 
\f— r wher^^r/ ~ {ms'^n/8C2)^^'^Ny/N. Dephasing of the condensate 
becomes significant only as A'^ — > A^, and not necessarily before the vortex 
lattice melts. 

However, in three dimensions, the Tkachenko mode frequency at long 
wavelengths becomes linear in the wavevector for any propagation direction 
out of the transverse plane.l^^ At zero temperature the vortex displacement 
correlations are convergent at large separation, but at finite temperatures, 
they grow with separation. The growth of the vortex displacements should 
lead to observable melting of vortex lattices at higher temperatures and 
somewhat lower particle number and faster rotation than in current exper- 
iments. At zero temperature a system of large extent in the axial direction 
maintains long range order-parameter correlations for large separation, but 
at finite temperatures the correlations decay with separation.'^ 
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6. BEYOND THE LLL REGIME 



At sufficiently liigh rotation, tlie vortex lattice should melt and become 
a vortex liquid. The regime just beyond melting has yet to be described in 
detail. At still higher rotation speeds, as seen in numerical simulations with a 
limited number of particles, the system begins to enter a sequence of highly 
correlated incompressible fractional quantum Hall-like states.'^'^^'^^'^ For 
example, at angular momentum Lz = 2N{N — 1), where A'^ (measured in 
terms of the total circulation, = {m/h)§v ■ di) equals 2N, the exact 
ground state is an A^-particle fully symmetric Laughlin wave function (in 
two dimensions), 

*(n,r2, . . . ,r^) ~ lliQ - 0)'e-^^'^i/24, (n) 

where Cj = Xj + iyj . Since the wave function vanishes whenever two parti- 
cles overlap, the interaction energy, ^gSj^j(5(rj — rj), vanishes in this state. 
Theoretically elucidating the states in general when the angular momentum 
per particle is of order the total particle number, as well as studying this 
regime experimentally, remain important challenges. 
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Fig. 3. Schematic phase diagram of the ground state of a rapidly rotating 
Bose condensed gas at zero temperature in the ri-Na/Z plane, where Z 
is the height in the z direction, showing the regions of multiply quantized 
vortices, of singly quantized vortices forming an array which at large number 
becomes a triangular lattice, and of an array of singly quantized vortices with 
a hole in the center. 



7. ANHARMONIC TRAPS 



The physics of a condensate confined in an anharmonic trap, e.g., Q, 
is quite different from that in a harmonic trap, since it becomes possible to 
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rotate the system arbitrarily fast. As the system rotates sufficiently rapidly, 
the centrifugal force pushes the particles towards the edge of the trap, and 
a hole open up in the center. Singly quantized vortex arrays with a hole 
have been seen in numerical simulations^^ and discussed theoretically in 
Refs. I1H32|33|34|351 In addition, at very high rotation, systems tend to 
form a single multiply quantized vortex at the center, with order parame- 
ter ~ e*'"^, where the integer quantization index u is >> 1. Such giant 
vortices have been seen in numerical simulations"^^ and are discussed the- 
oretically in Refs. I11I32I33I341 The schematic phase diagram, as a function 
of interparticle interaction strength vs. rotation rate is shown in Fig. 3. Full 
details can be found in Refs. I32I33I34I35I Initial studies of rapidly rotating 
condensates in harmonic lattices at the ENS are reported in Ref. IH71 
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